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ABSTRACT 

We investigate the formation of silicon carbide (SiC) grains in the framework of dust- 
driven wind around pulsating carbon-rich Asymptotic Giant Branch (C-rich AGB) 
stars in order to reveal not only the amount but also the size distribution. Two cases 
are considered for the nucleation process; one is the LTE case where the vibration 
temperature of SiC clusters is equal to the gas temperature as usual, and another is 
the non-LTE case in which is assumed to be the same as the temperature of small SiC 
grains. The results of hydrodynamical calculations for a model with stellar parameters 
of mass M*=1.0 Mq, luminosity L^,=10^ L©, effective temperature Tefj=2600 K, C/0 
ratio=1.4, and pulsation period P=650 days show the followings: In the LTE case, SiC 
grains condense in accelerated outflowing gas after the formation of carbon grains and 
the resulting averaged mass ratio of SiC to carbon grains of ~ 10^*^ is too small to 
reproduce the value of 0.01-0.3 inferred from the radiative transfer models. On the 
other hand, in the non-LTE case, the formation region of SiC grains is inner than 
and/or almost identical to that of carbon grains due to the so-called inverse greenhouse 
effect. The mass ratio of SiC to carbon grains averaged at the outer boundary ranges 
from 0.098 to 0.23 for the sticking probability as=0. 1-1.0. The size distributions with 
the peak at ~ 0.2-0.3 //m in radius cover the range of size derived from the analysis of 
presolar SiC grains. Thus the difference between temperatures of small cluster and gas 
plays a crucial role in the formation process of SiC grains around C-rich AGB stars, and 
this aspect should be explored for the formation process of dust grains in astrophysical 
environments. 

Subject headings: circumstellar matter; dust, extinction; stars: AGB and post-AGB; 
stars: winds, outflows 

1. INTRODUCTION 

Presolar silicon carbide (SiC) is the abundant and well studied dust species among presol 
ains extracted from meteorites. Based on the isotopic compositins of C, N, Si and Al, presol 
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SiC grains are classified into six categories named "mainstream", A, B, X, Y and Z grains (e.g., see 
Ottil201Cll for a review). Among them, the mainstream grains populating about 90 % of presolar 
SiC grains have been considered to be originated in carbon- rich Asymptotic G i ant Branch (C-rich 



AGB) stars from the isotopic compositions of C, N, and Si ( Zinner et al. 198'A Hoppe et al 



and the s-pr o cess s ignatures in the noble gases as well as in the heavy elements (jHoppe &: Ott 



19941 ) 



19971 ) 



Amari et al.l (jl994l ) have reported that SiC grains extracted from the Murchison meteorite have 
a log-normal distribution with the peak at ~ 0.4 //m in diam eter, and 70 % of th em by number 
are in the grains whose diameter ranges from 0.3 to 0.7 /im (jPaulton et al.l l2003l ) . Although it 
was reported that SiC grains are more fine grained in other meteorites than in Murchison (see 



Russell et al 



19971'). a sign ificant fraction of presolar SiC grains are larger than 0.3 ^m in diameter 



(see lHoppe &: Zinneijl2000l ) . Thus, submicron-sized SiC grains are expected to form around C-rich 
AGB stars. The presen ce of SiC graiii s in C-rich AGB stars was confirmed frorn the emission 
feature around 11.3 /um (lHackwelllll972l: iTreffers Cohenlll974l : iGoebel et al.lll98d ). prior to the 
discovery of presolar SiC grains (jBernatowicz et al.lll987l ). Al so, the absorption feature attributed 
to Si C grains was observed in several extreme carbon stars ()Jones et al.l Il978l : ISpeck et al.l 119971 . 
20091 ). The radiative transfer models fitting the observed spectral energy distributions (SEDs) 
including the spectral feature attributed to SiC grains have estimated the mass r a tio of SiC to 
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carbon gra i ns to be in the raiige of 0.01-0.3 (ILorenz-Martins &: Lefevrd 119931 . Il994l : iBlanco et al 



1199 



The variation of the mass ratio of SiC 
to carbon derived from the radiative transfer models arises from the grain models used, and the 
value could be also influenced by the optical constants of SiC used in the radiative transfer models 
(see the discussion in Section 6). 

Although the isotopic signatures of presolar SiC grains and the observed spectral feature around 
11.3 fxm support the formation of SiC grains in circumstellar envelopes of C-rich AGB stars, the 
formation process as well as the formation conditions remain unknown and debatable: Thermo- 
dynamic equilibrium calculations have been applied to iii vestigate the formation condition of dust 
grains in C-rich AGB stars (e.g. iLodders &: Feglevlll995l ). and the calculations have claimed that 
high gas pressure (0.3-300 dyne cm~^) and low C/0 ratio (1.0 < C/0 < 1.1) is necessary for 
reproducing the condensation sequence of presolar TiC, graphite, and SiC grains in this order in- 
ferred from the pre sence of graphitic sp herules containi ng a TiC core ([Sharp &: Wasserburgill995l : 
Daulton et al.ll2003l ). On the other hand. lMcCabd (jl982l ) suggested the formation of SiC grains ear- 
lier than carbon grains, considering the inverse greenhouse effect th at the temperature of SiC grains 



which is transparent in the optical to near infrared (NIR) region (jSpitzer 



195S 



Hofmeister et al 



20091 ) is lower than the temperature of carbon grains in circumstellar envelopes of C-rich AGB 
stars. In addition, the interferometric observation of IRC-l-10216 at 11 fiui has suggested the pres- 
ence of dust species more transparent th an carbon grains in the inner circumstellar envelope close 
to the photosphere (jPanchi et al.l Il990l ) . Although thermodynamic equilibrium calculations can 
tell us the physical conditions in which a condensate exists stably, the calculations as well as the 
astronomical observations provide directly no information on the number and size of grains that 
are essential to clarify the physical conditions prevailing in the formation sites of grains and the 
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physical and chemical processing suffering in the interstellar space. Also, it should be pointed out 
here that the C/0 ratio in C-rich AGB stars is not restricted to be less than 1.1, and tends to 
increase with decreasing the effective temperature Tgff; the low C/0 ratio observed in stars with Tes 
S, 2500 K is conside r ed to be caused by the neg lect of formation of dust in the model atmosphere 
(|Bergeat et ahlbood : Isergeat Chevanieil[2005f l. 



So far dust formation around C-rich AGB stars has been investigated in the fr amework of 
pulsation-enhanced dust-driven wind models since the study by [Fleischer et al.l (|l992l ). and it has 
been believed that the hydrodynamical models including only the formation of carbon grains well 
reproduce the obser ved dynamical behavior o f C-rich AGB stars w hose mass loss rates exceed 



10~^ Mq yr"^ (e.g.. IWinters et al.lll994 119971 : iNowotnv et al.ll2005l ). However, the formation of 



SiC grains has not yet pursued in the hydrodynamical models. On the other hand, formation of a 
variety of dust species including SiC grains has been studied in the steady state dust-driven wind 
models around AGB stars in order to investigat e the evolution of dust in the interstellar space (see 
Ferrarotti and Gailll2006l : IZhukovska et al.ll2008l ). However, in the models, the number of seed nuclei 
is assumed in the formation process. Although the models may predict the mass of dust species 
supplied from AGB stars, no information is available for the size distribution. Taking into account 
the nucleation process and including the so-called inverse greenhouse effect, iKozasa et al.l (jl996l ) 
have investigated the formation of SiC grains in an outflowing gas with constant velocity around 
C-rich AGB stars and shown that SiC grains condense prior to carbon grains for the mass-loss rate 
smaller than ~ 1.5 x 10^^ Mq yr~^; otherwise composite grains consisting of carbon and SiC form. 
The calculated radius of ~ 0.01 fim is too small to reproduce the size of presolar SiC grains and 
the proposed grains consist ing of a SiC core and a carbon mantle have not been discovered except 
for one in presolar grains (jCroat et al.ll201Cl ). Thus the investigation of formation of SiC grains 
in a more realistic hydrodynamical model remains an important subject to be explored which will 
reveal the formation process and condition around C-rich AGB stars. The knowledge of not only 
the amount but also the size distribution of dust gra. ins is crucial for investigating the role of dust 
in the universe as shown bv lYamasawa et al.l ((20111). In particular, the size distribution of SiC 
grains formed around C-rich AGB stars is vital to get insight into the evolution and processing 
during the journey from the formation sites to the incorporation into meteorites from the analysis 
of presolar SiC grains. 

In this paper, we aim at exploring the formation process of SiC grains and investigating 
the amount as well as the size distribution in the framework of pulsation-enhanced dust-driven 
wind from C-rich AGB stars. We develop the hydrodynamical model including for mation of car- 



bon a nd SiC grains as well as a non-grey radiative transfer. The formulation by iGauger et al 



(199C) is applied for the nucleation and growth processes of carbon grain^. We assert that SiC 



^In the formula tion it is assumed that carbon grains nucleate and grow homogeneously. The recent analysis of 
presolar graphites IjCroat et al.l 120051 ) reavealed that about 40 % of the carbide-containing graphites have a carbide 
core at their center. This implies that at least more than 60 % of presolar graphites nucleate and grow homogeneously. 
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grains nucleate and grow homogeneously starting from SiC molecule because all presolar SiC grains 
whose internal structure has thus far been analyzed do not contain any seed nuclei at their centers 



( Stroud &: Bernatowica |2005| ) . Furtheri nore there is no evidence in presolar graphite grains for 



a SiC seed nuclei except for one grain (jCroat et al.ll2010l ). The two cases are considered for the 



formation process of SiC grains; one is the LTE case where the vibration temperature of SiC 
cluster is the same as the gas temperature as usual. Another is the non-LTE case taking into 
account the inverse greenhouse effect on the formation, in which is assumed to be the same as 
the temperature of small SiC grains. We will show that the mass ratio of SiC to carbon grains and 
the radius of SiC grains are at most 10~^ and 0.1 iJ,m, respectively, in the LTE case, while the mass 
ratio ranges from 0.098 to 0.23 and the typical radius ranges from 0.2 to 0.6 ^m for the sticking 
probability Os of 0.1-1.0 in the non-LTE case. 

This paper is organized as follows; The chemical reactions and gas species relevant to the 
formation of SiC grains are examined and the formation process of SiC grains is formulated in 
Section 2, and the numerical scheme for calculation of grain formation is presented in Section 3, 
including the method for tracing the size distribution. Section 4 describes the hydro dynamical 
model and the modeling procedure for pulsation-enhanced dust-driven wind model including dust 
formation and radiative transfer. The results of calculations are presented in Section 5 and are 
discussed in Section 6. Summary is presented in Section 7. 



2. FORMATION PROCESS OF SIC GRAINS 



In the hydrodynamical models of pulsation-enhanced d ust-driven wind from C-rich AGB 
stars so far developed since the work of lFleischer et al.l (119921). the f ormat ion of carbon grains has 
been processed according to the recipe developed by ICauger et al.l (jl990l ): the abundances of gas 
species responsible for the nucleation and growth of carbon grains are assumed to be in chemical 
equilibrium in the gas w hose elemental abundance is solar except for carbon. In the recipe based on 
Gail &: Sedlmavii (|l988l ). the formation of clusters with number of constituent molecules (size) up 
to nout=1000 is evaluated from the steady-state nucleation rate, and the growth and evaporation 
of grains with size n > nout is treated as macroscopic grains. The formation of carbon grains are 
considered to proceed through the two body reactions of C2H2 and C2H molecules with carbon 
clusters, starting from atomic carbon or diatomic carbon molecule. 

In this section, we shall explore the formation process of SiC grains around C-rich AGB stars, 



follow ing the formulation for the formation of carbon grains (jCail &: Sedlmaviill988l : iGauger et al 
1990). First we examine the gas species responsible for the nucleation and the growth processes 



of SiC grains from chemical equilibrium calculations. In the same way as the formation of carbon 
grains, we consider that the nucleation process proceeds through the two body reactions of cluster 
with gas species, starting from the molecules whose stoichiometric composition is the same as SiC 
grains. On the other hand, we allow three or more body reactions on the surface for the growth 
process. Then, the basic equations describing the nucleation and growth processes of SiC grains 
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are presented. 



2.1. Si bearing gas species in dust formation region 



Figure 1 shows the temperature dependence of relative abundances of Si-bearing molecules 
(left) and C-bearing molecules (right) in chemical equilibrium in the gas with the density p = 
10"^'^ g cm~^ and C/0 = 1.4, which is calculated by the method of the min imizat i on of the Gibbs 
free energies. The elemental abund ances except f or car bon are taken fror nlAllenI (119731). and th e 



thermodynamic data are taken from lChase et al.l ( 19851 ) except for C2H2 ( Cherchneff et al 
SiCs, SiC4, SiCs, Si2C2, Si2C3, Si2C4, SiaC, Si3C2, SisCs, Si4C, Si4C2, and SisC (jPeng et al. 



199 



2), 



200. 



We can see that, as the starting molecules for the formation of SiC grains, SiC is most abundant 
in the region whose gas temperature is higher than 1350 K, while Si2C2 molecule is abundant in 
the lower temperature region. Among Si-bearing molecules that are considered to be candidate 
reactants in the growth process of SiC grains, Si, SiC2, and SiS molecules are abundant in the overall 
range of gas temperature presented in Figure 1, and Si2C is also abundant in the low temperature 
region of T < 1100 K. As shown in the right panel, hydrocarbon molecules such as C2H2 or C2H 
are expected to be reactants participating in the growth process on the surface of SiC grains unless 
these molecules are heavily depleted by the formation of carbon grains. 



2.2. Nucleation process of SiC grains 

Including the starting molecules SiC and Si2C2 for the nucleation process, the candidate re- 
actants for the two body reaction in the nucleation process are confined to twelve Si-bearing 
molecules: SiC, SiC2, SiCa, SiC4, Si2C, Si2C2, Si2C3, SisC, Si3C2, Si4C, Si4C2, and SisC. We shall 
exclude the gas species with which the reaction for formation of SiC clusters is endothermic. Table 

1 presents the reaction enthalpies /S.H° of the twelve candidates for the formation of Si2C2 and 
SisCs at 1000 and 1500 K. The molecules SiC2, SiCs, SiC4, Si2C, Si2C3, SisC, Si4C, and SisC are 
excluded from the candidates since the reactions are endothermic. The reaction enthalpies for the 
formation of SiC clusters with n > 3 for which no thermodynamic data is available are replaced 
with those for the formation of solid SiC, and the values at 1000 and 1500 K are listed in Table 

2 where the thermodynamic data of /3-Sicl^ is used in the calculations. The molecules SiC2 and 
SisC are excluded because of the endothermic reactions. Among the remaining candidates, as the 
reactants in the nucleation process of SiC grains, we adopt SiC and Si2C2 since these molecules 
are much more abundant than the others in the temperature region of interest to us as shown in 
Figure 1 (left). 



^ Presolar SiC grains so far analyzed show only two polytypic forms ( Daulton et al. 20031 ): /3-SiC (79.4 % by 
number), a-SiC (2.7 %), and the intergrowths of these two forms (17.1 %). Note that the difference of thermodynamic 
data between a-SiC and /3-SiC does not affect the results of calculations presented in this paper. 



- 6 - 



We discard the nucleation process starting from the two body reaction between Si2C2 molecules 
since Si2C2 is relatively stable and without the thermodynamic data of Si4C4 it is difficult to judge 
whether the reaction for the formation of Si4C4 proceeds or not. Thus we consider the following 
two reactions in the nucleation process of SiC grains: 



(SiC)„ + SiC 

(SiC)2„-l + (SiC)2 - 



> (SiC)„+i 

(SiC)2n+l 



(1) 

(2) 



where n > 1, and note that we take SiC as the starting molecule for the nucleation process of SiC 
grains. We shall discuss the effect of the nucleation process starting from the two body reaction of 
Si2C2 molecules in Section 6. 

In the nucleation process, clusters grow through the attachment of reactant and the growth 
rate depends on the number density and kinetic temperature of the reactants. On the other hand, 
clusters decay through the det achment of reactant ni olecules, and the de cay rate depends on the 
vibration states of the cluster (|Gail &: Sedlmaviill988l ). As pointed out bv lKozasa et al.l (|l996l ). in 
the dust formation region around C-rich AGB stars, the vibration temperature of SiC cluster 
is not always the same as the kinetic temperature of gas T due to the low gas density and strong 
radiation field. The principle of detailed balance being applied, the decay rate can be related with 
the growth rate by considering a hypothetical situation in which the cluster is in a local thermal 
equilibrium with the ambient gas whose temperature is T^. Then, under the assumption that is 
independent of the size of cluster, the steady-state nucleation rate Jg for the reaction with (SiC)j 
is given by 



Jl = 47rasaonsic?^? 



27rm,- 



1 + 



E 

^=1 {ik + i; 



2/3 



(3) 



where i=l (2) denotes the reaction (1) (the reaction (2)) and the quantities related with the 
reactant SiC (Si2C2) molecule, oq is the hypothetical radius of a monomer in the bulk phase, 
nsic the number density of starting molecule SiC. The number density, sticking probability, mass, 
and partial pressure of the i~th reactant is represented by n^, a*, rrii, and Pi, respectively. The 
Boltzmann constant is ks, -/Vx = [(jT-out ~ 2)/iJ with the floor function [x\, and the surface tension 
of the cluster of size n an- The supersaturation ratio Si{T^, Pi) is given by 

ln5,(r.,P,) = -^+ln#. (4) 



where Pq is the standard pressure. In the above equation, AG^ is the Gibbs free energy for 
the formation of bulk SiC from the i-th reactant in the standard pressure, and the temperature 
dependence of AG-'/Zcb^v is evaluated from the least-squares fitting given by 

AG? 



ksT^ 



ai + ^ilnT^ + 



Vi 

+ — + — H — -■ 



(5) 
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whos e form is based on a polynomial expression of heat capacity at the constant pressure (jBerman &: Brown 
19851 ). The coefficients Oj, 7,, 5i, rji and ej are given in Table 3. In the calculations, we adopt 
ao=l. 71x10"^ cm which is evaluated from the bulk density of solid SiC (/3=3.16 g cm~^) taken 
from the CRC hand book (jLidd l2003l ) . The values of the sticking probabilities of SiC and Si2C2 
molecules on the SiC cluster are unknown and are assumed to be unity. The surface tension (T„ of 
SiC cluster with size n=2 or 3 is directly evaluated from the following definition of surface tension 
of n-mer 



47rao(ra 



1)3CJ„ 



g{n) -°g{l) - {n - l)g^ 



(6) 



where g{n) and Qs are the Gibbs free energies of n-mer in gas phase and the monomer in bulk 
phase in the stand ard state, respectively. O n the other hand, we adopt the bulk value of C7=840 erg 
cm~^ measured by lAllen Kingeryl ( 1953 ) as the surface tension of SiC clusters with size n > 3. 



2.3. Growth, of SiC grains 



We allow three or more body reactions for the growth of SiC grains in addition to the two 
body reactions for the nucleation. As shown in Figure 1, the abundant Si~ and C~bearing molecules 
expected to be the reactants for the grain growth are Si, SiC2, Si2C, C, C2H2, and C2H, except 
for CO and SiS which are stable against the chemical reactions. Then we select the following 
exothermic reactions; 



(SiC)„ + SiC ^ (SiC)„+i 
(SiC)„ + Si2C2 ^ (SiC) n+2 
(SiC)„ + 2Si + C2H2 ^ (SiC)„+2 + H2 
(SiC)„ + 2Si + C2H + H ^ (SiC)„+2 + H2 
(SiC)„ + 2Si2C + C2H2 ^ (SiC)„+4 + H2 
(SiC)„+Si + SiC2 ^ (SiC)„+2 
(SiC)„+Si + C^ (SiC)„ 



+1 



(7) 
(8) 
(9) 
(10) 

(11) 
(12) 

(13) 



where n > nout=1000; the reaction enthalpies at 1000 K are -1265, -1739, -1572, -1206, and -1240 
kJ mol~^ for the chemical reactions from (9) to (13) in this order. 

Except for the two body reactions, it is assumed that each of the chemical reactions pro- 
ceeds through the attachment of a key species which is defined as the species of the least colli- 
sion frequencies among the reactants and is consider ed to control the kinetics of growth process 
(jKozasa &: Hasegawalll987l : iHasegawa k. Kozasalll988l ). Then, the time evolution of the size of dust 
grains with radius a is given by 



d 
dt 



1 



3rn 



(14) 
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where the size-independent net time scale of grain growth Tnet (jGauger et al.lll99d ) is the sum of 



the net time scales for the reactions from (7) to (13) and is written as follows: 

In the above equation, r^g^ (i=l,2) is the net time scale of grain growth by the two body reaction 
with (SiC)j molecule and given by 



1 ; 2 i kBT f _fev \ 

^— = z^7rar,a^\ rii ; , (16) 

and r^'" (n=9 to 13) is the net time scale of grain growth for each of the chemical reactions from 
(9) to (13) and is given in the same form of equation (16) by replacing a*, rrii and rii with those 
for the key species. In equation (16), Pev represents the vapor pressure for the key species in the 
reaction; for example, in the reactions (9) and (10), the vapor pressure of Si atom is assumed to 
be the same as that for the reaction: (SiC)s — )• Si + C. In the reaction (11), the vapor pressure 
of Si2C molecule is assumed to be the same as that for the reaction: 2(SiC)s — s- Si2C + C. In the 
calculations, the sticking probability is assumed to be unity for all reactions. 



3. NUMERICAL SCHEME FOR GRAIN FORMATION 

The number densities of the gas species responsible for the formation process of carbon and 
SiC grains in the pulsation-enhanced dust-driven wind from C-rich AGB stars are evaluated from 
the chemical equilibrium calculations for given gas density and temperature coupled with the hy- 
drodynamical calculation. We include 4 atoms (H, C, Si, S) and 14 molecules (H2, C2, C3, C2H, 
C2H2, SiH, SH, II2S, SiC, SiC2, Si2C, Si2C2, CS, SiS) in the chemical equ ilibriu r n calc ulations, and 
the abundances of elements except for carbon are taken from the table of lAllenI (|l973l ). We employ 
th e formulation present ed in Section 2 for the formation process of SiC grains, while the formulation 
bv ICauger et al.l (1990j) is applied for the formation of carbon grains. As mentioned before, the 
sticking probability ag for all reactions at nucleation and growth of SiC grains is assumed to be 
unity in the calculations, whose effect on the results of calculations is discussed in Section 6. 



The number density and size of grains as well as the fraction of condensible molecules locked 
into grains are evaluated by taking the j'-th moment of the size distribution of grains per H-element 



Kj defined as 



n=nout 



(17) 



where f{n,t) is the number density of grains with the size n at time t and n<-H>=^H+ 2nHT Nor- 

mally, the moments are calculated by reducing to the simultaneous differential equations (jCail &: Sedlmavi 
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19881 ). The time variation of gas (grain) temperature is not always monotonic and rather compli- 
cated in the pulsation-enhanced dust-driven winds. In the case that the formed grains are being 
lost by the destruction due to evaporation at time t, the extinction term is included in the differ- 
ential equations by using the number density of grains with the minimum size /(j^outji) which has 
to be derived by tracing the size distribution. Therefore, solving the differential equations must be 
compensated with the evaluation of the size distributio n of grains with a limited number of size 
bins in the hydro dynamical simulations (jFleischeii Il994l ) . 



Thus, in this paper, we calculate the moments by employing the finite and adaptive size grids 
rather than so lving a set of the siinultaneous differential equ ations. Each of the size bins is comoved 
in size space (jKriiger et al.l Il995l : IWoitke &: Niccolinil |2005| ) . The i-th size bin is characterized by 
the number of grains per H-element w{i) in the size interval with lower limit 

di^rnin Sind Upper limit 

Oj.max- The lower and upper limits are calculated by integrating equation (14). When the net time 
scale Tjiet is positive, the number of grains per H-element in the lowest size bin 'w{I) where / is the 
total number of size bins is given by integrating the equation 



dw{I) Js 



dt n<H> 

When the net time scale Tnet is negative, with the destruction rate of grains Jjes estimated by 

da t()(/)n<H> 



(18) 



Jde 



dt fl/, upper Q-lowest 



(19) 



1/3 

where aiowest = ^o^-outi '"^(-^) is derived by integrating the equation 

dw{I) Jdes 



dt 



?^<H> 



(20) 



According to lWoitke Sz Niccolinil ( 20051 ). the discrete moment Kf^^'^^ is evaluated by 



'7-1 



K 



discr 



(j + l)«o 



E 

A:=0 



owest'^/.max 



(21) 



k=0 



for j=0 to 3. The 0-th moment K^^^'^^ and the 3rd moment K^^^'^^ give the number of grains per 
H-element and the number of condensible molecules incorporated into grains per H-element, respec- 
tively. The volume equivalent radius of grains (a) at a time is defined by (a) = ao(i^f ''7i^o''")^^^• 



4. THE MODEL 



4.1. Hydrodynamics 



The basic equations for the time-dependent spherically symmetric hydrodynamical model are 
written in standard finite difference form and solved by explicit integration along the lines of 
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difference scheme give n bvlRiclitmver &: Mortonl (1196711. This modeling method has been adopted 



by B erhn group (e.g., [Fleischer et al.lll992l : [Winters et al.l l200d : iHelling et al.l l200d : iDrever et al 



201 ll ) . Actually, this method is time-consuming since the tiine step s are restricted by the Courant 



Friedrichs-Levy (CFL) stability condition (see iHofner et al.lll996l ). However, since the advection 



term drops out from the equation of motion in Lagrangian coordinate, we can easily trace the size 
distribution of grains which is one of main subjects of this study. 

The continuity equation representing the conservation of mass is replaced with the definition 
of gas velocity v given by 

dr 

where r is the radial position of the Lagrangian mass element at a time t. Under the assumption 
of position coupling between gas and grains, the equation of motion is described as 

dv _ I dp GM{r) 
dt p dr 



(22) 



a 



(23) 



where p is the thermal gas pressure, p is the gas density, G is the gravitational constant, and M{r) 
is the mass inside a radius r. In the above equation, a is the ratio of radiation pressure force to 
gravity and expressed as 



a 



47rcGM(r 



(24) 



:as 
pr 



where Kpr is the sum of the flux weighted mean of mass radiation pressure coefficient for gas K 
and that for dust Kp"*^*, L^, the luminosity, and c the speed of light. The mass inside a radius 
r in the above equations is replaced with the stellar mass because the mass contained in the 
circumstellar envelope is small compared to the stellar mass during the calculations pr esented in th is 
paper. In the calculations, we adopt the constant gas opacity of «;g=2xl0~^ cm^ g~^ (jBowenlll988l ). 



Following the previous works for hydrody namic models of dust-driven win d solved with frequenc y 



dependent radiative transfer calculation (iHofner et al. 



(Hofner et al. 


2003: 


Woitke 


2006: 



Mattsson et all l2010l ) 



we evaluate the dust opacity k!^^^ at a wavelength A in the Rayleigh approxi mation. The o ptical 
constants of carbon a. nd SiC grains used in the calculation are taken from Draine (|l985l n and 
Choyke &: PalikI (jl985l ). respectively. The radiation field in the circumstellar envelope is derived 
from the frequency-dependent radiative transfer calculation using the variable Eddington factor 
method under the assumption of radiative equilibrium. Then, the gas temperature Tgas is given by 



,gas 
^A,abs 



B\{T^as)d\ 



(25) 



■^The Rosseland m ean and Planck mean absorption coefficients calculated by the optical constants of 
aphite ( Drainel 1985 ) are almost the same as those calculated by the optic al constants of amorphous carbon 



^Rouleau fc Martini 1991 ) used in the recent hydrodynamical calculations (e.g., Hofner et al. 20031 : Mattsson et al 
I2OIO). Thus, the difference in the optical constants does not so much affect the results of calculations presented in 
this paper. 
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where the subscript "abs" refers to absorption, J\ is the mean intensity and B\{Tgf^s) the Planck 
function. The temperature of i-th grain species Tdust,j is given by 

/ '^tab^^'^AC^A = / t^t,lhsB\{Tdnst,i)d\. (26) 
JO JO 



4.2. Vibration temperature of SiC cluster 



Accord ing to the hydrodynamical models devel oped so far for dust-driven wind from C-rich 
AGB stars (jFleischer et al.lll992l : iHofner et al.lll995l ). the vibration temperature of carbon cluster 
is set to be equal to the gas temperature as usual. On the other hand, in order to investigate the 
influence of the inverse greenhouse effect suggested by iMcCabd (jl982l ') for the formation of SiC 
grains, we consider the two cases for the vibration temperature of SiC cluster: One is the case that 
the vibration temperature of SiC cluster is the same as the gas temperature and that is referred 
to as the LTE case hereafter. Another is the non-LTE case whe re the vibration tem perature is 
assumed to be the same as the temperature of small SiC grains (jKozasa et al.lll996l ). Assuming 
that the gas species mainly consisting of H, H2, and He accommodate with grains completely upon 
colliding, and placing a small SiC grain hypothetically at a position in the gas flow, we derive the 
vibration temperature of SiC cluster by balancing the heating with the cooling through the 
interaction with radiation and g follows; 



A!abs(«cl)[-^A 



"-gas J 




gas 



T 



(27) 



where QA'abs("ci) is the absorption efficiency factor of SiC grain with radius ad whose value is set 
to be 10~^ /im in the calculation, and ngasj, rrij^ and Cyj are the number density, the mass, and 
the specific heat at constant volume of the j-th gas species, respectively. 



4.3. Modeling procedure 

The models of pulsation-enhanced dust-driven wind are specified by the six parameters; the 
stellar mass M^,, the stellar luminosity L*, the effective temperature Teff, the C/0 ratio, the period 
of the pulsation P, and the velocity amplitude of the pulsation Aup. As the first step of the 
hydrodynamical simulation, a dust-free static model atmosphere is constructed as follows; The 
initial stellar radius i?* is given by the Stefan-Boltzmann law = AirB^aT^^ where a is the 
Stefan-Boltzmann constant. Given a guess value for gas density at i?*, the radial profiles of 
density and temperature are calculated by solving the equation 

1* = (.8) 

p ar 
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coupled with the radiative transfer calculation with sufficiently fine grids. The procedure is iterated 
until the geometrically diluted optical depth at i?* converges to 2/3. 

The dynamical model calculations including the formation of grains are carried out by placing 
the piston at i?ino whose position is two times pressure scale height around i?* below from the 
photosphere of static model atmosphere and is defined as the inner boundary of dynamical model 
at t=0. Then the radius Rf^^^ of inner boundary is varied sinusoidally with time as 

i^fr^ = Rino + Aup£ sin (l^t) (29) 

and the velocity uf^^'^ at the inner boundary is given by 

= AnpCos(^t). (30) 

In the simulation, no mass flow is assumed through the inner boundary. The radiative flux is 
assumed to be constant over time and be equal to 'itB{Tcs) = vr Bx{Tci[)dX at the photosphere 
of dynamical model -Rpht where the diluted optical depth of gas is 2/3, and thus the stellar luminosity 
L=K = Aj^Rp^^aT^g varies with time. The piston velocity amplitude Aup is increased slowly from 1 
cm s~^ to the specified value to prevent t he first outwardly moving subsonic wave from growing 



into a enormous transient shock, following iBowenI (jl988l ). 



After the onset of dust formation, the dust-driven wind takes places. As the gas expands, the 
spacing between radial grids in the acceleration region becomes wider and the rezoning procedure 
necessary for resolving the density structures of the circumstellar envelope is introduced by adopting 
the method proposed in Fleischer et al.(1992). After the outermost zone passes through the outer 
boundary of dynamical model which is placed at 25 i?* in the simulation, the outermost Lagrangian 
grid is eliminated. Then, at the outer boundary, we trace the time variations of mass-loss rate, gas 
outflowing velocity, dust-to-gas mass ratio, and condensation efficiency fc (/si) which is defined as 
the fraction of carbon (silicon) atoms in C-bearing molecules except for CO (Si-bearing molecules) 
locked into carbon (SiC) grains. Note that, in the case that all Si atoms in Si-bearing molecules 
are locked into SiC grains, the maximum value of fc ~ 0.9 in the model with C/0=1.4. 



5. RESULTS 



We adopt the following model parameters in the calculation; Tcff=2600 K, L*=10^ Lq, M=f=1.0 
Mq, C/0=1.4, P= 650 days, and A^p=2.0 km s'K whi ch is a typical set of the values in the pre- 
vious studies (e.g., iFleischeii Il994l : IWinters et al.l l200d ) and is employed as the reference values 
i n the models of C-rich AGB stars that reason ably reproduce the observed dynamical behaviors 
(jWinters et al.lll994l . 119971 : iNowotnv et al.ll2005l ). The number of radial grids used in the hydrody- 
namical calculations depends on the density structure and is in the range between 800 and 1200 
throughout the calculations. The wavelength region covered in the radiative transfer calculation 
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ranges from 0.35 to 120 fxm and the number of wavelength grids is 100. The simulations have been 
performed for 180 cycles (stellar pulsation periods). In this paper, we neglect the radiation pressure 
force acting on SiC grains because of the transparency in the optical to infrared region (up to ~ 9 



/im) (|Chovke Paliklll985l : iPitman et al.1 l2008l : iHofmeister et al.ll2009l ) 



The formation of SiC grains in the pulsation-enhanced dust-driven wind is closely related with 
the formation of carbon dust shell stemming from the rapid formation of carbon grains in the high 
density gas induced by the stellar pulsation and the resulting acceleration of dense region caused 
by the radiation pressure force acting on the carbon grains. The carbon dust shell is characterized 
by the inversions of density distributions of gas and carbon grains and plays a crucial role in the 
time evolution of gas temperature and density through the backwarming caused by the thermal 
emission from carbon grains. First, we outline the formation of carbon dust shells in this model, 
and then we show how SiC grains form in the LTE and non-LTE cases. 



5.1. Formation of carbon dust shells 

The carbon dust shell appears in the present model quasi-periodically once a stellar pulsation 
period and is classified into two types by the peak value of fc in the shell. Here we shall overview 
the formation of each type of the carbon dust shells and its dynamics, referring Figure 2 which 
displays the radial structure of the quantities related to gas dynamics and formation of carbon and 
SiC grains in the LTE case throughout the period from t=166.7 to 168.8 P. 

The first type of carbon dust shell (hereafter referred to as the 1st type CDS) is defined as 
the high density region with the peak value of fc ^ 0.9, and is located around 4 R^: at t=167.9 P 
(see the bottom row in Figure 2). The formation of this 1st type CDS actually starts at t ~ 166.7 
P and the behavior of formation process is almost the same as that of the 1st type CDS whose 
formation starts at t ~ 168.8 P; the nucleation of carbon grains takes place in the dust-free region 
ranging from 1.8 to 3.0 i?* whose gas temperature decreases to 1300 K (see the radial variation 
of the nucleation rate of carbon grains at t=168.8 P and the time variation of the distribution 
of the number of carbon grains in this region from f=168.5 to 168.8 P). Then, the backwarming 
by thermal emission from carbon grains subsequently formed in the outer part of this region not 
only raises up the gas temperature but also evaporates the carbon grains in the inner part (see the 
3rd rows at t=167.0 P and f=167.3 P). However, in the outer part, the nucleation and growth 
of carbon grains still proceed and the ratio of radiation pressure force to gravity a around 2.5 i?* 
exceeds unity at t=167.0 P with increasing the volume equivalent radius (ac) and the condensation 
efficiency fc of carbon grains up to ~ 0.1 ^um and ~ 0.3 (see the 3rd and bottom rows), as a result 
this part is accelerated outwards. On the other hand, the gas density in the region outer than 
the accelerated region is insuffcient to produce abundant enough carbon grains to accelerate this 
region efficiently despite the ongoing nucleation. Thus, the radial compression due to the velocity 
difference causes the outer part of the accelerated region to accelerate more through the increase 
of nucleation and growth rates of carbon grains, and results in the shell structure; note that the 
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shell structure is not associated with the pulsation shock but induced by the formation of carbon 
grains. The shell structure becomes prominent with decreasing the gas density in the region inner 
than the accelerated region due to the radial expansion (see the 2nd rows from t=167.3 to 167.6 
P). The formation of 1st type CDS completes at t ~ 167.9 P. 

The formation of the second type of carbon dust shell (hereafter referred to as the 2nd type 
CDS) defined as the high density region with the peak value of 0.3 < /c ^ 0.7 onsets around 2.2 
i?* at t=168.2 P, being triggered by the merging of the two consecutive shocks induced by the 
stellar pulsations, which can be seen from the time variations of gas velocity and a in this region 
from t=167.9 to 168.2 P; the carbon grains in front of the precedent shock, which nucleated at i ~ 
166.7 P and survived from the evaporation process associated with the evolution of 1st type CDS 
starting at t ~ 167.0 P, grow up from (ac)=0.05 to 0.1 fim in the region compressed further by 
the merging with the following shock (see the bottom rows from t=167.6 to 168.2 P). After then, 
the nucleation and growth of carbon grains in the compressed outer part of this region develops 
the shell structure with increasing a, while the rapid decrease of gas density caused by the outward 
velocity already exceeding 10 km s^^ suppresses the growth rate of carbon grains and the increase 
of fc- The evolution of this carbon dust shell ends up as the 2nd type CDS with the peak value 
of fc ~ 0.7 at t ~ 169.0 P. The density profile as well as the peak value of fc of 2nd type CDS 
varies from cycle to cycle (e.g., see the 2nd type CDS located around 4.5 at t=167.0 P). The 
alternating formation of the 1st and 2nd type CDSs is sometimes disturbed, in which case the 2nd 
type CDS fails to grow up without increasing the peak value of fc up to 0.3. The appearance 
and periodicity of carbon dust shell de p end on the stellar p aram eters in a complica t ed in anner 



as investigated by Fleischer et al.l (119951 ). iHofner et al.l (119951 ) . and lDreyer et al.l (|2009l . 1201 ll ). For 



example, in our model with C/0=1.8, only the 1st type CDS forms every stellar pulsation period. 



5.2. Formation of SiC grains in LTE case 

The values of mass-loss rate, gas terminal velocity, and dust-to-gas mass ratio averaged over 
the last sixty pulsation periods at the outer boundary are 1.24x10"^ Mq yr-i, 20.4 km s'^, and 
1.64 xlO^'^, respectively, in the hydrodynamical model with the LTE case for the formation of SiC 
grains. As can be seen from Figure 2, the nucleation of SiC grains onsets when the gas cools down to 
1100 K after the formation of carbon grains, and is activated in the temperature range of 800-1000 
K. The corresponding formation region is located around 3.5 ~ 5.5 R^, unless the gas temperature 
is raised up higher than 1100 K by the backwarming from the precedent 1st or 2nd type CDS. 
In the formation region, Si2C2 molecule is more abundant than SiC molecule and the nucleation 
proceeds through the reaction (2) rather than the reaction (1). In particular, the formation of SiC 
grains is most efficient in the region around the outer boundary of the 1st type CDS because of 
the lower temperature and higher gas density, and the nucleation rate per H-element reaches the 
maximum value of ~ 10~^^ s~^ at i=167.9 P. However the subsequent decrease of the gas density 
caused by the accelerated radial expansion depresses the nucleation and growth processes of SiC 
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grains; the number of SiC grains per H-element is limited to smaller than 10~^^ (see the bottom 
row at t=168.8 P), and the volume equivalent radius of SiC grains (asic) is typically 3-6x10"^ 
fim. In the other regions, the number and size of SiC grains are much smaller. Thus the depletion 
of Si-bearing molecules does not take place substantially in LTE case (see the radial variation of 
/si in the bottom row throughout the time), and almost all of SiC grains form around the outer 
boundary of carbon dust shell. 

Figure 3 represents the the radial distribution of gas and grains in the region of r < 25 R^, 
at t=166.8 P. The radial distribution of fc (the 2nd row) clearly shows that the different type of 
carbon dust shell forms by turns as addressed in Section 5.1. The carbon dust shells broaden with 
moving outwards from the formation sites because the outer part with more abundant carbon grains 
is accelerated more efficiently than the inner part, and the broadening in the outer circumstellar 
envelope of r > 10 i?* is caused by the pressure gradient. In addition, the slower moving 2nd type 
CDS seems to be caught up by the faster moving 1st type CDS around 21 R^. Note that the shell 
(shock) structure seen in the radial distribution of gas density precedes the corresponding carbon 
dust shell in the outer circumstellar envelope (see the top row). The distributions of psic (the 
top row) and /si (the 2nd row) demonstrate that SiC grains reside in thin shell around the outer 
boundary of carbon dust shell, reflecting the formation site. The value of fsi (psicZ/Ogas) and the 
mass ratio of SiC to carbon grains SiC/C are at most 10~^ (10~^) and 10~^, respectively (see the 
lower three rows). The value of fc {pc/Pgas) averaged over the last sixty pulsation periods at the 
outer boundary is 0.707 (1.64 x 10^^), while the averaged value of fsi (/Osic//Ogas) is in the order of 
magnitude of 10~^ (10~^^). The resulting averaged value of SiC/C is 4.93 x 10~^ which is much 
smaller than the value inferred from the radiative transfer models (0.01-0.3). 

Figure 4 shows the size distribution of SiC grains by mass averaged over the last sixty pulsation 
periods at the outer boundary. The distribution is characterized by the single-peaked profile with 
the peak around a=6.0xl0~^ pm, and the radius of SiC grains is limited to be less than 0.1 pm. 
The modal size is almost a fa ctor of four small er than the size of presolar SiC grains extracted 



from the Murchison meteorite (lAmari et al .11 19941 ). The LTE case with lower C/0 ratio resulting in 



lower gas outfiow velocity makes the formation of SiC grains with radius of 0.2-0.3 pm possible as 
long as the gas outfiow exhibits a layered structure with the CDSs. On the other hand, regardless 
of C/0 ratio, the amount of SiC grains formed in the LTE case is too small to reproduce amount 
of SiC grains inferred from astronomical observations. 



5.3. Formation of SiC grains in non LTE case 

Figure 5 displays the radial structure of the quantities related to gas dynamics and formation 
of grains in the non-LTE case throughout the period from t=166.7 to 168.8 P. The vibration 
temperature of SiC cluster defined in Section 4.2 and depicted in the top row (solid red) is evaluated 
by placing the small SiC grain with aci=10^'^ pm hypothetically in the dynamical calculations. The 
difference between the gas temperature Tgas and the vibration temperature is negligible in the 
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high density region of p > 10~^^ g cm~^ where the colhsion with gas dominates the energy balance 
and controls T^. As the gas density decreases and the collision with gas gets less effective than 
the interaction with the radiation in the energy balance, deviates from and turns to be lower 
than Tgas be cause small SiC grain i s almost transparent from the visible to in frared region (up 



to ~ 9 fim) (jChovke k PalikI Il985l : IPitman et al.1 boosi : iHofmeister et alj|2009l ). Also, it should 



be noted that is not affected so much by the backwarming due to the thermal emission from 
the outer carbon dust shell compared with the gas temperature. Thus, the temperature difference 
AT = Tgas — 7v is magnified in the less dense region inner than the opaque carbon dust shell (e.g., 
see around 2.5 at t=167.3 P). 

As the gas density decreases, the supersaturation ratio of SiC grains quickly exceeds unity with 
increasing AT as well as decreasing T^. The nucleation process of SiC grains can be activated in 
the region with p < 10~^^ g cm~^ as long as Tv < 1200 K and AT > 200 K, even in the region with 
Tgas ^ 1400 K where the nucleation process of carbon grains is almost depressed. On the other 
hand, if is higher than 1200 K, the nucleation process of SiC grains is not activated regardless 
of the gas temperature since the abundances of molecules responsible for formation of SiC grains 
are not sufficient to make the supersaturation ratio large enough for the onset of nucleation. 

The active nucleation and growth of SiC grains takes place in the infalling gas around 1.8-2.2 
i?* at t=166.7 P when the gas density drops down to a few times 10"^'^ g cm~^, Tv decreases down 
below 1200 K and AT increases by larger than ~ 200 K , almost concurrently with the onset of 
formation of the 1st type CDS as explained in Section 5.1. Then, the backwarming from the evolving 
1st type CDS raises up Tv and Tgas high enough so as to evaporate SiC grains in the inner dense 
region with p > 10^^'^ g cm~'^ at t=167.0 P. However, in the outer less dense region, SiC grains 
grow up to ~ 0.3 ^m thanks to the increase of AT and the condensation efficiency /si reaches 0.1. 
In the region around 2.5 i?*, the nucleation of SiC grains proceeds, and at t=167.3 P the increase 
of fsi is distinguishable while carbon grains evaporate due to the backwarming, which continues up 
to t ^167.5 P (see the 3rd and bottom rows). On the other hand, in the region around 2 R^,, the 
evaporation of SiC grains turns to the growth at t=167.3 P with weakening the backwarming as the 
1st CDS moves outwards, and then the increase of gas density up to 10^^^ g cm^^ caused by the 
pulsation shock leads to the efficient nucleation and growth of SiC grains, and /si increases up to ~ 
0.7 quickly until t=167.6 P. Then, the efficient growth of SiC grains in the region compressed by 
the merging of consecutive shocks starting from t ~ 167.9 P results in the complete consumption 
of Si-bearing molecules (/si=1.0) around t ~168.1 P, and (asic) in this region is typically larger 
than ~ 0.2 pm, in contrast to the LTE-case. The region with the condensation efficiency fsi ~ 
1.0 is referred to as the Si-depletion region hereafter. Concurrently, the 2nd type CDS formed and 
developed in this region stretches the Si-depletion region outwardly (see the radial distribution of 
/si in the bottom row from i=168.2 to 168.8 P), while the backwarming from the 2nd type CDS 
makes T^ higher than 1200 K and depresses the nucleation of SiC grains in the region inner than 
the Si-depletion region for a while (see Tv in the top row at i=168.5 P). It should be noted that 
the 2nd type CDS moves outward more slowly than the 1st type CDS as mentioned in Section 5.2. 
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Thus, the Si-depletion region with slower outward velocity is finally caught up and pushed by the 
1st type CDS as can be seen from the time evolution of the Si-depletion region spreading over from 
3.5 to 5.5 iZ* at t=167.3 P which formed two cycles before. 

As mentioned above, the formation of Si-depletion region onsets in the infalling gas in front 
of the coming pulsation shock, and the region moves outward with the formation of 2nd type CDS 
triggered by the merger with the next pulsation shock. Thus, the Si-depletion region forms once 
every two pulsation periods. Also, even in the other regions, the resulting value of /si often exceeds 
0.1. Thus, in contrast to the very small value of at most 10~^ in LTE case, a large amount of SiC 
grains can condense through the inverse greenhouse effect in the non-LTE case. 

Figure 6 presents the radial variations of gas and grains at t=176.3 P. It seems that SiC 
grains are concentrated around the outer boundaries of 1st type CDSs in the density distribution 

(sec the top row) because the shock associated with the formation of carbon dust shell sweeps up 
the gas. Reflecting the variation of gas density, the radial distribution of psjc docs not always 
corresponds to the Si-depletion region shown in the 2nd row. It should be noted that the extended 
Si-depletion region does not develop without the 2nd type CDS, which can be demonstrated from 
an example of the failure of alternating formation of both types of carbon dust shells located at 
18 i?* where Si-depletion region restricted in a narrow region resides just ahead of the 1st type 
CDS formed about 1 cycle after from the formation of the Si-depletion region (see the 2nd row). 
Apart from the 1st type CDS located at 18 i?*, the extended SiC dust shell corresponding to the Si- 
depletion region between 1st and 2nd type CDSs can be clearly recognized in the radial distribution 
of Psic/Pgas (see the 3rd row where /9sic/Pgas=8.46xlO~^ corresponding to /si = 1.0). The radial 
distribution of /si shows that the Si-depletion region stretched over by the slower moving 2nd type 
CDS narrows with moving outwardly, being compressed by the faster moving 1st type CDS. Thus, 
in the outer circumstellar envelope, the radial width tends to be of the order of the stellar radius 
(sec /si around 24 i?* where the 2nd type CDS almost merges with the 1st type CDS) and the 
narrower Si-depletion region seems to reside around the outer boundary of the 1st type CDS. 

The radial distribution of the mass ratio SiC/C has a sharp peak in the region just ahead 
of the outer boundary of 1st type CDS, and decrease by 0.1 (10~^) toward the outer boundary 
of precedent 2nd (1st) type CDS (see the bottom row). The peak value exceeds ten in the inner 
circumstellar envelope, and decreases with increasing r because of the subsequent growth of carbon 
grains. However the value of SiC/C in the Si-depletion region still keeps larger than 0.2 with the 
peak value exceeding unity. 

The resulting values of mass-loss rate, terminal velocity, dust-to-gas mass ratio averaged 
over the last sixty stellar pulsation periods are 1.10x10"^ Mq yr-\ 19.5 km s-\ and 1.85x10"^, 
respectively, for the non-LTE case. The reduction of terminal velocity compared to the LTE case 
reflects the somewhat weakened acceleration resulting from the reduction of amount of carbon grains 
due to the incorporation of C-bearing gas species into SiC grains; on average 4.40 % of carbon 
atoms in the condensible C-bearing molecules is incorporated into the SiC grains. The averaged 
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value of pc/Pgas (/c) is 1.51 x 10"^ (0.652), while the value of psic/pgas (/si) is 3.41 xlO"'' (0.403). 
The resulting averaged value of SiC/C is 0.226 and is close to the upper value derived from the 
radiative transfer models. The size distribution of SiC grains by mass averaged over the last sixty 
stellar pulsation periods at the outer boundary is depicted in Figure 7. The size distribution shows 
the broad profile with the mean radius of 0.3 /ira, in contrast to the sharp and single-peaked size 
distribution in the LTE case presented in Figure 4. In the size range from 0.03 to 4.0 /xm, 80 % 
of SiC grains by mass populate in the radii between 0.2 and 1.0 /xm: The distribution well covers 
the size range of presolar SiC grains found in the Murchison meteorite. Thus the hydro dynamical 
model with the non-LTE case for the formation of SiC grains that the vibration temperature of SiC 
cluster is assumed to be the same as the temperature of small SiC grains reasonably reproduces the 
amount of SiC grains inferred from astronomical observations and the size of presolar SiC grains 
extracted from the Murchison meteorite. 



6. DISCUSSION 

In the calculations of hydrodynamical models presented in Section 5, we assume that SiC is 

the starting molecule for the nucleation process of SiC grains, excluding the process starting from 
the two-body reaction of Si2C2 molecules. This may underestimate the amount of SiC grains in 
the LTE case since Si2C2 molecule is more abundant in the lower temperature region than SiC 
molecule as shown in Figure 1 (left). In the LTE case where carbon grains condense in prior to SiC 
grains, the nucleation process starting from the two-body reaction of Si2C2 molecules dominates 
the nucleation rate in the region around 6-8 i?* whose gas temperature ranges from 600 to 800 
K. However the gas density in the region of p ~10~^^ g cm~^ is too low to increase the amount 
of SiC grains; even if this process is included in the calculation, the resulting values of both /gi 
and SiC/C are limited to 10~^ although the nucleation rate is much larger than that starting from 
SiC molecule. In addition, this nucleation process produces the population of grains whose radii 
are much smaller than 6.0 x 10~^ /xm, and make the modal radius smaller than 3.0 x 10~^ fim. 
Thus the nucleation process starting from two-body reaction of Si2C2 molecules does not change 
the conclusion that the hydrodynamical models with the LTE-case cannot reproduce the amount 
of SiC grains around C-rich AGB stars inferred from the astronomical observation and the size 
range derived from the analysis of presolar grains. 

In the non-LTE case, this nucleation process dominates the formation of SiC grains in the 
region of r >4 R^, and is activated particularly in the outer part of the carbon dust shell. Although 
the number of SiC grains produced through this process reaches about 10 times that produced 
through the nucleation starting from SiC molecule, the typical radius of SiC grains formed in the 
region with low gas density is limited to 3x 10~^-6x 10~^/im. Thus, the inclusion of this nucleation 
process does not affect the size distribution of SiC grains whose radius is larger than 0.1 fim. Also 
it is true for the amount of SiC grains since the averaged /si increases by only a few percent. 
Thus we conclude that the exclusion of nucleation process starting from the two-body reaction of 
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Si2C2 molecules does not affect the result of calculations in the non-LTE case, except for the size 
distribution in the range of asic < 6xl0~^/um. 

The sticking probability set to be unity for nucleation and growth processes of SiC grains 
may overestimate the size and amount of SiC grains formed in the non-LTE case. Figure 8 shows 
the time averaged size distributions by mass of SiC grains for given values of ag. In the range of 
as=0. 1-1.0, the mass fraction of large sized grains with a > 1.0 //m decreases significantly and 
the profile becomes narrow with decreasing Os, while the reduction of does not affect the mass 
fraction of small-sized grains with a <0.1 /xm. The modal radii are 0.30 and 0.22 fim for as=0.5 
and 0.1, respectively. Although the SiC grains with radii larger than 1.0 //m do not form for the 
case of as=0.1, the hydrodynamical models in the non-LTE case rea sonably reproduce the size 
range of presolar SiC grains extracted from the Murchison meteorite (jAmari et al.lll994l ) as long 
as Us > 0.1. The averaged value of SiC/C decreases with decreasing as; the values are 0.199 and 
9.77x10"^ for 0^=0-5 and 0.1, respectively, which is within the range of the value inferred from 
the radiative transfer models. Although the reduction of changes the radial distribution of SiC 
grains substantially, we conclude that the hydrodynamical models in the non-LTE case cover the 
ranges of size derived from the analysis of presolar grains and amount evaluated from the radiative 
transfer models unless the sticking probability as ^0.1. 

The mass ratio of SiC/C inferred from the radiative transfer models depends on the models 
used in the calculation; the evaluat ed values for C-rich AGB stars w ith the emi ssion feature of SiC 
grains are in the range of 0.01-0.29 ( Lorenz-Martins Lefevr3ll994 ). 0.01-0.15 ( Groenewegen 



Groeneweeen et al 



199 



and 0.1-0.3 (Blanco et al 



199 



The smaller value by iGroenewegen 



1995: 



(I1995I ) and IGroenewegen et al.l (|l998l ) is considered to be ascribed to the grain model; the tem- 



perature of composite grain consisting of carbon and SiC used in their mode l is hi gher than that 
of isolated SiC grains used in the other models as is shown in iBlanco et al.l (|l998l ). since carbon 
grains are more absorptive in the visible to NIR region than SiC grains. Therefore, in order to 
reproduce the observed SEDs, if the inner boundary of dust layer is placed close to the photosphere, 
the high temperature of com posite grains causes the va l ue of SiC / C to reduce. Also, the smaller 
value of 0.01-0.06 derived by iLorenz-Martins &: Lefevrd (119941 ) for C-rich AGB stars showing SiC 
emission feature with optically thick shells (ri^m > 3.0) would arise from the optical constant of 
SiC grains used in the models; the radiative transfer models usually employ the optical constant s 
derived from the SiC powder synthesized in the laboratory (iBorghesi et al.lll985l : lPegourig|l988l ). 
The absorption coefficients in visible t o infrared reg i on (up to ^ 9 iim) is much larger than the 
crystalline SiC (jChovke &: PalikI Il985l : iPitman et al.l 120081 : iHofmeister et al.l |2009| ) , which results 
in the higher temperature of SiC grains and the smaller SiC/C. In addition, the mid-infrared ob- 
servation of C-rich AGB stars by Spitzer shows that the 30 /xm emission feature attributed to 
MgS grains becomes prominent with decreasing the strength of SiC e m ission feature in the star s 
with mass-loss rate larger than 10~^ Mq yr~^ ( Leisenring et al.l 120081 ). IZhukovska &: Gaill ( 20081 ) 
demonstrates that the absorption efficiency around 11.3 of grains consisting of a SiC core and 
a MgS mantle decreases with increasing the volume fraction of MgS mantle. Therefore, if it is true, 
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the radiative transfer model could underestimate the value of SiC/C without including SiC grains 

coated by MgS mantle. Thus, the mass ratio of SiC to carbon grains derived from the the radiative 
transfer models depends on the optical constants of SiC as well as the model of grains used in the 
calculations, and the value of SiC/C up to ~ 0.23 calculated in the non-LTE case for C-rich AGB 
stars is considered to be not always in conflict with the values derived from SED fittings. 

7. SUMMARY 

The formation of SiC grains around C-rich AGB stars is investigated for the first time in the 
framework of hydrodynamical model for the pulsation-enhanced dust-driven wind. We formulate 
the nucleation and growth processes of SiC grains, considering that SiC grains nucleate and grow 
homogeneously starting from SiC molecule. In the calculations, the two cases arc considered for 
the nucleation process of SiC grains; the LTE case in which the vibration temperature of SiC 
cluster is equal to the gas temperature and the non-LTE case in which Ty is assumed to be the 
same as the temperature of small SiC grain whose radius is 10~^ /xm. On the other hand, the 
vibration temperature of carbon cluster is considered to be the same as the gas temperature as is 
usually assumed in the hydrodynamical models for C-rich AGB stars. 

The results of calculations for the model parameters M*=1.0 M©, L*=10^ L©, reff=2600 K, 
C/0 ratio=1.4, P=650 days, and Aitp=2.0 km s~^ are summarized as follows: In the LTE case, SiC 
grains form in the accelerated and low-density outflowing gas after the formation of carbon grains. 
The resulting time averaged number fraction of Si atoms locked into SiC grains (condensation 
efficiency /si) of at most 10~^ is too small to reproduce the amount of SiC grains inferred from the 
astronomical observations and the radius is limited to less than 0.1 /xm. On the other hand, in the 
non-LTE case, the time averaged mass ratio of SiC to carbon grains ranges from 0.098 to 0.23 for 
the sticking probability as=0. 1-1.0, which is not in conflict with the value of 0.01-0.3 inferred from 
the radiative transfer models. The time averaged size distribution of SiC grains by mass with the 
peak at radius of 0.2-0.3 /xm well covers the size range of presolar SiC grains extracted from the 
Murchison meteorite. Thus we conclude that the so called inverse greenhouse effect plays a crucial 
role in the formation process of SiC grains in the pulsation-enhanced dust-driven winds from C- 
rich AGB stars, apart from the validity of the underlying assumption made in the non-LTE case. 
The non-LTE effect should be explored for the formation process of dust grains in astrophysical 
environments. 

The hydrodynamical model provides the two types of carbon dust shells; the 1st (2nd) type 
CDS with the peak value of /c ^0.9 (0.3 ^ /c ^ 0.7). The formation region and the resulting radial 
distribution of SiC grains in the hydrodynamical model with non-LTE case is closely related with 
the formation and dynamics of carbon dust shell, especially through the density enhancement and 
the backwarming effect. In the inner circumstellar envelope, the SiC dust shell can be discriminated 
from the 1st type CDS through the extended Si-depletion region caused by the 2nd CDS formed 
in the outer part of active SiC formation region. However, in the outer circumstellar envelope of 
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r >20 i?*, the SiC dust shell is localized around the outer boundary of 1st type CDS. Not only the 
formation region and the radial distribution but also the amount and the size distribution of SiC 
grains could heavily depend on the C/0 ratio as well as the other model parameters; for instance, 
in the non-LTE case with M*=2.0 Mq, the mass ratio of SiC to gas and the size of SiC are almost 
the same as those with M*=1.0 Mq while the mass loss rate decreases by a factor 5. Anyway 
more comprehensive study covering a wide range of model parameters is necessary for revealing 
the amount and the size distribution of SiC grains formed around C-rich AGB stars. Also, the 
difference between the radial distributions of carbon and SiC grains could greatly influence the 
appearance of the spectral feature attributed to SiC grains. These aspects will be investigated in 
the forthcoming papers. 
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Fig. 1. — Relative abundances of Si-bearing molecules (left) and C~bearing molecules (right) in 
chemical equilibrium in the gas with the density p = 10"^'^ g cm~^ and C/0 ratio = 1.4. 
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Fig. 2. — Radial structure in LTE case for i=166.7 to 167.6 P (top panel) and for t=167.9 to 168.8 

P (bottom panel). Top row; the gas velocity v (blue) and the gas temperature Tgas (black). 2nd 

row; the gas density p (blue) and the ratio of radiation pressure force to gravity a (black). 3rd 

row; the nucleation rates per H-element Js/^<H> of carbon (black solid) and SiC (red) grains, 

the destruction rate per H-element Jdes/''^<H> of carbon grains (black dashed), and the numbers 

per H-element Kq^^" of carbon (grey) and SiC (green) grains. Bottom row; the condensation 

efficiencies fc (black) and /si (red), and the volume equivalent radii (a) of carbon (grey) and SiC 

(green) grains. 
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t=166.8P 




■/R* 



Fig. 3. — Radial distributions of gas and grains at t=166.8 P in the LTE case. Top row; the 
mass densities of gas (blue, Pgas), carbon (black, pc) and SiC (red, psic) grains. 2nd row; the 
condensation efficiencies /c (black) and /si (red). 3rd row; the density ratios of pc/Pgas (black) 
and psic/Pgas (red). Bottom row; the mass ratio of SiC to carbon grains SiC/C. 
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Fig. 4. — Size distributions of SiC grains by mass averaged over the last sixty stellar pulsations at 
the outer boundary in the LTE case. Note that g{a) is the size distribution function and the value 
of g{a)a^ is normalized to the maximum value. 
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t = 166.7 P t = 167.0 P t = 167.3 P t = 167.6 P 
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Fig. 5. — Same as Figure 2, but for the non-LTE case. Note that the red hue in the top rows 

denotes the vibration temperature and the red-dashed hne around 1.9 i?* in the 3rd row at 

t=167.0 P depicts the destruction rate of SiC grains per H-element Jdes/"'<H>- 
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t=176.3P 




r/R* 

Fig. 6. — Same as Figure 3 but for t=176.3P in the non-LTE case. 
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Fig. 7. — Same as Figure 4, but for the non-LTE case. 
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Fig. 8. — Averaged size distributions of SiC grains by mass for given values of as, each of which is 
normaUzed by the maximum value. 
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Table 1: Reaction enthalpies Ai7° for growth reactions of SiC cluster with size n=2 to 3. 







AH° [kJ mor^] 


Reaction 


lOOOK 


1500K 


SiC + SiC ^ SiaCa 


-751.5 


-751.0 


SiC + SiaCa - 


> SiaCg 


-518.8 


-514.5 


SiCa + SiC ^ 


Si2C2 + C 


76.98 


78.12 


SiC2 + Si2C2 - 


SigCs + C 


309.6 


314.7 


SiCs + SiC ^ 


Si2C2 + C2 


-54.88 


-58.76 


SiC3 + Si2C2 - 


Si3C3 + C2 


177.8 


177.8 


SiC4 + SiC ^ 


Si2C2 + C3 


-123.0 


-131.1 


SiC4 + SisCs - 


->■ Si3C3 + C3 


109.7 


105.4 


Si2C + SiC ^ 


Si2C2 + Si 


-110.8 


-109.8 


Si2C + Si2C2 - 


Si3C3 + Si 


121.9 


126.7 


Si2C2 + Si2C2 


^ Si3C3 + SiC 


232.7 


236.5 


Si2C3 + SiC - 


> Si2C2 + SiC2 


-199.5 


-202.6 


Si2C3 + SiC - 


> SiaCg + C 


110.1 


112.1 


Si2C3 + Si2C2 


^ Si3C3 + SiC2 


33.14 


33.97 


SisC + SiC -> 


Si2C2 + Si2 


-53.73 


-52.46 


SisC + Si2C2 - 


Si3C3 + Si2 


178.9 


184.1 


Si3C2 + SiC - 


> Si3C3 + Si 


-122.4 


-119.9 


Si4C + SiC -> 


Si2C2 + Si3 


-122.6 


-124.9 


Si4C + Si2C2 - 


Si3C3 + Si3 


110.1 


111.6 


Si4C2 + SiC - 


> Si2C2 + SigC 


-195.5 


-199.5 


Si4C2 + SiC - 


> Si3C3 + Si2 


-16.53 


-15.41 


Si4C2 + Si2C2 


^ Si3C3 + SigC 


37.19 


37.04 


SisC + SiC ^ 


Si2C2 + Si + Si3 


352.2 


347.9 


SisC + Si2C2 - 


-^Si3C3 + Si+Si3 


584.8 


584.5 
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Table 2: Reaction enthalpies Ai7° for growth reactions of solid SiC. 







AH° [kJ mol-i] 


Reaction 




loUUK 


SiC + SiC(s) — 


► 2SiC(s) 


-793.8 


-788.5 


S1C2 + SiC(s) - 


2biC(s) + C 


34.64 


40.66 


SiCs + SiC(s) - 


2SiC(s) + C2 


-y I.2Z 


-yo. zz 


SiC4 + SiC(s) - 


2SiC(s) + C3 


-165.3 


-168.6 


SisC + SiC(s) - 


-> 2SiC(s) + Si 


-153.1 


-147.3 


Si2C2 + SiC(s) 


^ 3SiC(s) 


-836.2 


-826.0 


Si2C2 + SiC(s) 


^ 2SiC(s) + SiC 


-42.34 


-37.47 


Si2C3 + SiC(s) 


^ 3SiC(s) + C 


-207.2 


-199.4 


SigCg + SiC(s) 


^ 2SiC(s) + SiC2 


-241.9 


-240.0 


SigC + SiC(s) - 


2SiC(s) + Si2 


-96.07 


-89.92 


Si3C2 + SiC(s) 


^ 2SiC(s) + Si2C 


-286.6 


-284.1 


Si3C2 + SiC(s) 


^ 3SiC(s) + Si 


-439.8 


-431.3 


Si4C + SiC(s) - 


2SiC(s) + Sia 


-164.9 


-162.4 


Si4C2 + SiC(s) 


^ 2SiC(s) + SiaC 


-237.8 


-237.0 


Si4C2 + SiC(s) 


^ 3SiC(s) + Si2 


-333.9 


-326.9 


SisC + SiC(s) - 


2SiC(s) + Si + Si3 


309.8 


310.4 



Notes. SiC(s) denotes solid SiC which is assmued to be /3-SiC. 



Table 3: Coefficients for calculating AG°/A;bTv in equation (5). 



i 


ai 


A 


7i 


5^ 


Vi 


1 


3.11656x10^ 


-1.12418x10° 


1.71784x102 


-1.01082x10^ 


1.01432x10^ -8.78548x10^ 


2 


-2.96847x10^ 


2.79219x10^ 


5.55156x10"^ 


-8.91964x10^ 


1.15759x10*^ -1.28618x10^° 



